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Consider two normal leads coupled to a superconductor; the first lead is biased while the second one and 
the superconductor are grounded. In general, a finite current /2(Vi,0) is induced in the grounded lead 2; its 
magnitude depends on the competition between processes of Andreev and normal quasiparticle transmission 
from the lead 1 to the lead 2. It is known that in the tunneling limit, when normal leads are weakly coupled 
to the superconductor, Ja(Vi) 0) = 0, if j Vi < A and the system is in the clean limit. In other words, Andreev 
and normal tunneling processes compensate each-other. We consider the general case: the voltages are below 
the gap, the system is either dirty or clean. It is shown that 72(Vi,0) = for general configuration of the 
normal leads; if the first lead injects spin polarized current then I2 = 0, but spin current in the lead-2 is finite. 
XISIN structure, where X is a source of the spin polarized current could be applied as a filter separating spin 
current from charge current. We do an analytical progress calculating Ii(Vi, V2), ia(Vi, X^). 

PACS: 74.50.+r, 74.80.Fp, 75.75.+a 



Hybrid systems consisting of a superconductor (S) 
and two or more normal metal (N) or ferromagnetic 
(F) probes recently started to attract a great atten- 
tion ^EJIjUQJE]. Among most striking new results is 
the prediction that NSN (FSF) devices can play the role 
of entangler producing Einstein-Podosky-Rosen (EPR) 
pairs 0] having potential applications, for example, in 
quantum cryptography 6 . Not long ago rather un- 
usual effect was described in normal metal - tunnel 
barrier (I) - superconductor - tunnel barrier - nor- 
mal metal (NISIN) junction (see, e.g., Fig. lb) 0[3]. 
It was shown that when Ni is biased, N2 and S are 
grounded there is no current injection from Ni to N2 
at subgap biases; main assumptions were: 1) the su- 
perconductor is clean, 2) large number of conducting 
channels are involved in electron tunneling through NS 
interfaces [21131 ■ In other words, the subgap cross con- 
ductance G12 = dv 1 l2{Vi,0)\\v 1 \<A = 0, where the cur- 
rent Ji flows in Ni, Vi is the bias between Ni and S 
and V2 — between N2 and S. The suppression of G12 
was attributed to the compensation of the contributions 
to the current from Andreev and normal quasiparticle 
tunneling processes between Ni and N2 PJ- It was also 
noted that G12 7^ in FISIF junctions: G12 decays 
exponentially as exp(— r/£) with the characteristic dis- 
tance r between the normal terminals (see, e.g., Fig^i), 
where £ is the superconductor coherence length; at small 
r /£,, G12 decays also rather quickly (at atom-scales): as 
l/(fcF?") 2 (&f m the superconductor) [2|. Thus with 
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Figl. The outline of the setup. Ni,2 are normal metals 
or ferromagnets. 



clean superconductors a measurement of G12 may be- 
come difficult. 

In this letter we first of all generalize results [3131 
and get rid of the assumption 1) [i.e., S is not restricted 
to be clean] . We show that when the superconductor is 
dirty (the mean free path is smaller than £) Andreev and 
normal transmission rates [as well as G12 in FISIF junc- 
tions] slowly decay with the characteristics distance r 
between the normal (ferromagnetic) terminals (at r < £) 
in contrast to the clean regime (see Refs. [HE])- ^or 
example, in FISIF with superconducting layer thinner 
than £, see Fig. \iy>, G12 ~ ln(r/£); when the supercon- 
ductor is bulk then G12 ~ £/r [r > \p is supposed]. 
Measurements of the effects, related to electron tunnel- 
ing through a superconductor (e.g., G12) in the dirty su- 
perconductor case can be easier realized experimentally 
then in the clean case because then r is not restricted to 
atomic scales but by £ 3> Xf- We show that contribu- 
tions to the current from Andreev and normal quasipar- 
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Fig2. Electron scattering from a NSN junction. 



ticle tunneling processes always compensate each other 
in NISIN junctions (so, e.g., I 2 (Vi,V 2 = 0) = for 
| V\ | < A in first nonvanishing order over the transparen- 
cies of the layers I) for any amount of disorder in the 
S-layer. If one prepares a NISIN junction with layers I 
having large transparency then normal tunneling start 
dominating Andreev tunneling (and 12(^1,^2 = 0) 7^ 
0). We also considered FISIN junction, in particular 
with Vf 7^ and Vjv = 0. Then the ferromagnet F plays 
the role of the spin-polarized current injector. In this 
case I2{Vf, Vn = 0) = also, but spin current in N is fi- 
nite: charge component of the current converts into the 
supercurrent, spin accumulates in N. So XISIN struc- 
ture, where X is a source of the spin polarized current, 
could be applied in spintronics as a filter separating 
spin current from charge current. We find Andreev The 
and normal transmission probabilities T ee of a NISIN 
sandwich for subgap energies \E\ < A and different an- 
gles 9 between incident quasiparticle trajectory and the 
normal to NS interface. It is shown that the probabil- 
ities have resonances where The ~ T ee ', averages of The 
and T ee over incident channels (over 9) are equal — this 
is the reason why ^2(^1, V<i = 0) is suppressed and the 
spin current l!f\Vi, V% — 0) is finite. 

We start investigation of NISIN structures from the 
sandwich sketched in Fig^,: barriers at NS boundaries 
provide spectacular reflection; electrons in N and S move 
ballistically; the number of channels at both NS bound- 
aries is much larger than unity. The transmission prob- 
abilities T he (E, 9) and T ee {E, 9) [see Fig|2] describe An- 
dreev and normal tunneling of an electron incident on 
the NS boundary with the angle 9 and the energy E 
correspondingly into a hole and an electron in the lead 
2. Following the Landauer-Buttiker approach |gll9HlT)|: 

h{V U V 2 = 0) = ~J dE J2 ( T ee ~ The)x 

channels 

x(/ (1) -/ (2) ), (1) 

where the sum is taken over channels (spin degrees 
of freedom are included into channel definition); f^ 1,2 ' 
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Fig3. Resonances of Andreev and normal transmis- 
sion probabilities T ee (9), Th e (6) of a NISIN junction 
Fig^i at different transparencies Tjvs(#) of the layers 
I. Parameters: d/£ = 0.1, E F /A = 1000, the energy 
E — 0. Resonances correspond to k F dcos(8) /n — n, 
n = 1,2, .... If < E < A then the shape of the 
graphs slightly change: the resonance peaks become 
slightly asymmetric [same applies to the case when Tns 
are different but small]. It can be checked (even ana- 
lytically) that the areas under corresponding resonance 
peaks of T ee and The become equal at small Tns- 



are distribution functions in the leads 1,2; e.g., — 
riF^E) = 1/[1 + exp(/3E)], f( l > is not necessary a Fermi- 
function. We calculate the transmission and reflection 
probabilities using Boguliubov equations (BdG). The 
layers I are approximated by (5-barriers. Quasiparticle 
motion parallel and perpendicular to the NS interfaces 
can be decoupled |11II12| . Matching appropriate wave 
functions in the normal region and the supercondctor 
we get 8x8 linear system of equations for the transmis- 
sion amplitudes. Analytical progress can be made. It 
follows that if there is no barrier at NS boundaries (ex- 
cept A) The/Tee < (A/£ F ) 2 for any thickness d of the 
superconducting layer. This result is intuitively quite 
clear because A <C Ep can hardly reverse the direction 
of the quasiparticle momentum being about Uf [TT1u3| . 
However if there are barriers at NS boundaries in ad- 
dition to A (e.g., insulating layers I) then the situation 
changes: at certain 9 the transmission probabilities have 
resonances where The ~ T ee . When the transmission 
probabilities of the layers I, T N $ <§C 1, the areas un- 
der the resonance peaks of T e h(9) and T ee (9) are nearly 
same and 



(The) ~ (?ee), 



(2) 
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where (...) = E c „(---)/^„ * Jo(- ■ O^os 2 9. 
Eq. J5J is exact in first nonvanishing order over Tns- 
The resonances appear at kFdcos(9 n ) — nn, n = 
1,2..., give the leading contributions to (The), (Tee) 
and are responsible for Eq.J2J. The resonance width 
T ~ min{l, T/vs, d/£}. Typical dependencies of T e h(9) 
and T ee (9) from 9 and T N $ are illustrated in Fig|2J 

In fact 9 is discrete variable; its particular value 
is determined by the channel of the incident parti- 
cle. Equation is applicable when 1) Tns{9) slightly 
change when 9 changes from one channel to an adjacent 
one and 2) change of 9 from one channel to another 
should be smaller than the resonance width. The con- 
dition 1) is fulfilled typically when Tns(9) -C 1, the 
condition 2) requires Af/VA <C min{l, Tns(0)> d/£}, 
where A is the junction surface area. 

It follows from Eqs. illlL'li that subgap charge injec- 
tion from the lead 1 into the lead 2 in weak coupling 
regime (Tns <C 1) is suppressed: h(Vi,V2 = Q) = 0, 
because charge currents of transmitted hole and elec- 
tron quasiparticles compensate each other in the lead 
2; all the electron current converts into Cooper-pair su- 
percurrent in S. However if spin-polarized current is in- 
jected from the lead 1 finite spin current appears in 
the lead 2; transmitted electron and hole quasiparti- 
cles contribute the spin current. XISIN structure with 
Tns ^ -Tivl ^ 1 (this condition allows to neglect the 
contribution to the charge current going in S from An- 
dreev reflection at NiS surface), where X is current "in- 
jector" , can play the role of the filter of spin and charge 
currents, see fig.0] Equation for the spin current follows 
from Eq. JTJ: 

4 s) (Vi,V 2 = 0) = ^ fdE <7i{T ee + T he )x 

channels 

x(/ (1) -/ (2) ), (3) 

where a\ — ±1 labels spin degrees of freedom in X. 
General feature of transmission probabilities T and the 
current — their exponential suppression with d/£ when 
d ^> £ (£ is the superconductor coherence length). 

We show below that all the results discussed above 
remain true in general NISIN structure with more com- 
plicated shape than in Fig|TJi, (e.g. like in Fig. QJd) no 
matter dirty or clean. 

In general a system of weakly coupled normal (ferro- 
magnetic) and superconducting layers can be described 
by the Hamiltonian: H = Hi + H2 + Hs + Ht, where 
H1.2 refer to the electrodes Ni and N 2 , and Hs to the 
superconductor. The tunnel Hamiltonian Ht, which 
we consider as a perturbation, is given by two terms 
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Fig4. Isolator - superconductor - isolator- normal 
metal structure can help to separate in space spin and 
charge components of the current. 



Ht = H^ + H}p' corresponding to one-particle tunnel- 
ing through each tunnel junction: 
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where indices i = 1, 2 refer to normal (ferro) electrodes; 
tjj.p is the matrix element for tunneling from the state 
k = (k, a) in normal lead Nj to the state p = (p, a 1 ) 

a (i) 

in the superconductor. The operators d k and b p cor- 
respond to quasiparticles in the leads and in the super- 
conductor, respectively. 

The current can be expressed through the quasi- 
particle scattering probabilities within the Landauer- 
Biittiker approach. It is possible to calculate the scat- 
tering probabilities within the tight-binding model J3J, 
but it is more convenient to describe the current on 
the language of electrons only: Andreev transmission 
probability The Q is closely related to the crossed An- 
dreev (CA) tunneling rate r^[(Vi, V2) which shows how 
many electron pairs tunnel per second from the leads 1 
and 2 into the condensate of the superconductor (each 
lead gives one electron into a pair) and vice versa cor- 
respondingly, see Fig. [Sp, and [2]. Elastic co-tunneling 
rate ^ec) corresponds to T ee . Direct Andreev tunnel- 
ing rates, an< ^ ^(DA) describe Andreev reflec- 
tion in the leads 1 and 2 [see, e.g., Fig.[S^]. The current 



in the lead 2 consists of two contributions: one, L 



2 ' 



comes from the electron injection from the lead 1 due to 
crossed- Andreev and cotunneling processes, the other, 
I2 , ~ from the direct electron-tunneling between the 
lead and the superconductor. Same applies for the lead 



1. Thus h{V u V 2 ) = I^iVx, V 2 , -r i 2 



4 I '(v 1 ,v 2 ) = r {EC) -r (CA) , 



l2 D '(V 2 ), where 



T (EC) 
T (CA) 



4 D) (v 2 ) 



1 (EC) ~ 

1 (CA) 

r S^2 
1 (DA) ~ 



1 (EC)' 

L (CA) ' 

r S^2 
(DA) ■ 



(5a) 
(5b) 
(5c) 
(5d) 



Using the Fermi Golden rule, the rates can be found 
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Fig5. a) Direct Andreev tunneling (Andreev reflection), 
b) Crossed Andreev tunneling (Andreev transmission), 
and c) Elastic Cotunneling (normal transmission). 



in second order in the tunneling amplitude tk,p. Follow- 
ing the approach described in Ref. [21 ll4lfT5| . we finally 
obtain 



■E^(2VA^e), (6) 



[a 2 - e 

where ny) is the distribution function in the lead i — 
1,2. Hereafter we take H = 1, e = 1 [we do not assume 
ny) to be only equilibrium Fermi function]. The rate 
^CA 12 can b e obtained from the expression for r^7 12 
by substitution of (1 — n) for n. 

The kernel E { f A) (s) = dt~£ A) {t)e~ st is the 

Laplace transform of E^ A \t). It can be expressed 
through the classical probability P(Xi,pi;X 2 ,p 2 ,t) 
meaning that an electron with the momentum directed 
along pi initially located at the point X\ near the NS 
boundary arrives at time t at some point X 2 near the NS 



boundary with the momentum directed along p 2 spread- 
ing in the superconducting region as follows 

4 CA) (*) = 

= 87r 3g4^ s J dpi,2 J dX lj2 P(X 1 ,p 1 ;X 2 ,p 2 ,t)x 

(6{X u X 2 y 



G^(X 1> p 1) a)G ( - 2 HX 2) p 2) a) sin 2 f 
+ G (1) {X 1 , pt , g)G {2) {X 2 , pa, -<r) cos 2 ^ 



0(*i,* 2 ) 



Here the spatial integration is performed over the NiS 
and N 2 S surfaces. We choose the spin quantization axis 
in the direction of the local magnetization in the termi- 
nal N 1 (2)- The quasiclassical probabilities G^(X,p,a), 
i = 1,2 for the electron with spin polarization a tun- 
neling from the terminal N< to the superconductor are 
normalized in such way that the junction normal con- 
ductance per unit area g& (X) and the total normal 
conductance G^f are determined as [T4l[l3) 



= dpG^(X,p,a), G 



W 

N 



Then the normal conductance per unit area, discussed 
above, is defined as g$ = G$ /A, where A is the sur- 
face area of the junction. Symbol 9(Xi,X 2 ) is the angle 
between the magnetizations of the terminals Nj and N 2 
at points X\ and X 2 near the junction surface. If elec- 
trons in Ni and N 2 are not polarized then 9 = 0. 
In a similar way: 

= ^ 3 / dexx 1 ^ - w - 4 2 >(e - 

A 2 



[a 2 - e 2 j 

where 

4^)(t) = 

= j^3~4^ / d Pi^ 2 / dX 1 , 2 P(X 1 ,p 1 ;X 2 ,p 2 ,t)x 
{G^\X u p u o)G^{X 2 ,p 2 ,-a)^ 
+ G« (A"i,pi, cr)G^ (X 2 ,p 2 , a) cos 2 (^^2 ) | 



(7) 



The rate ^e~c) can be obtained from the expression 

(EC) 



for by substitution of (1 — n) for n. DA-rates 



are written in |14|. Equations (|5al7|l derived here al- 
low to describe transport properties of many types of 
junctions. 
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Consider FISIN junction with biased ferromagnet 
with the respect to the superconductor, the normal 
metal N has same voltage as S. So the ferromagnet 
plays the role of a current "injector" ; electrons com- 
ing from F are distributed with some distribution func- 
tion n^. Electrons in the deep of the terminal N 
are Fermi-distributed. It follows from Eas. (|5al7|l that 
contributions to the current from EC and CA pro- 
cesses compensate each other for subgap voltages so 
In(Vf 7^ 0,0) = 0. However spin current is finite: 

I^(V F ,0)=A^ f d( :J2<r[nP(t-V F )-nW(0] 

X [J^-e]BJe*v 3 J dPl ' 2 J dXl * X 
P(X 1 ,p 1 ;X 2 ,p 2 ,t)G^(X 1 ,p 1 ,a)G^(X 2 ,p 2 ). (8) 

Finally we consider a FISIF junction. It was 
shown in that in this junction I 2 (Vi,0) ^ and 
I 2 (Vi , 0) changes its sign when the ferromagnetic termi- 
nals change their orientation from parallel to antiparal- 
lel. Naively it can be supposed that in a FISIN junction 
where F is a current injector, S, N are grounded spin ac- 
cumulation at the interface of the normal metal would 
lead to spin-splitting of the density of states in N and a 
charge current. However this is not so, this corrections 
are of higher order over tunneling amplitudes than the 
processes in Fig. and can be neglected because we 
assume that tunneling amplitudes are small. 

It was also noted in that the cross conductance 
G12 = 9^/2(^1,0)1^=0 is suppressed in a FISIF struc- 
ture as l/(kpr) 2 when the characteristic distance be- 
tween the ferromagnets r < £. In dirty regime there 
is no conductance suppression at atomic-scales. Con- 
sider, for instance, the layout sketched in Fig. QJj; the 
width d of the superconducting film is supposed to be 
smaller than £. According to Eas. l|5al7|) the cross- 
conductance dependence from the distance r is deter- 
mined by the Laplace transform P(s = 2yA) of the 
probability P(r,t) = cxp(-r 2 / 4D\t\) /4nd\t\, where D is 
diffusion constant in the superconductor, d < £. When 
A F < r < f, G12 ~ P ~ hi(r/£) and if r > f, 
G12 ~ P ~ exp(— r/£). When the superconductor is 
bulk (d > £) similarly we find G\ 2 ~ %/r, Af <C f < £• 
All considerations above apply also for CA- and EC- 
rates. Thus it is practically more convenient to mea- 
sure finite effects related to electron subgap tunneling 
through a superconductor when it is dirty rather than 
clean. In dirty case the terminals are not restricted to be 
as close as Xf like in clean case but closer then £ ^> Xp- 
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